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SUMMARY 


Measurements  of  cloud  particle  images  and  concentrations  using 
laser  spectrometers  housed  in  pods  mounted  on  the  wings  or  fuselages  of 
research  aircraft  can  be  affected  by  the  distortion  of  the  airflow  about  the 
aircraft.  The  Flight  Research  Laboratory  has  developed  a  pod  with  a  Rose- 
mount  858  5-hole  probe  and  pressure  transducers  to  measure  airflow  angles 
and  velocities  at  typical  mounting  locations  on  cloud  physics  research 
aircraft.  This  report  documents  the  results  of  an  extensive  wind  tunnel 
calibration  of  this  pod  to  determine  the  factors  relating  the  differential 
pressure  measurements  to  the  flow  angles  and  velocities,  and  in  particular 
to  account  for  the  effects  of  the  canister  itself  on  these  measurements. 


RESUME 


Les  mesures  d’images  et  de  concentrations  de  particules  de  nuage 
au  moyen  de  spectometres  a  laser  loges  dans  des  nacelles  montees  sur  les 
ailes  ou  le  fuselage  d’aeronefs  de  recherche  peuvent  etre  biaisees  par  des 
perturbations  du  mouvement  de  l’air  autour  de  Faeronef.  Le  Laboratoire 
de  recherche  en  vol  a  mis  au  point  une  nacelle  munie  d’une  sonde  Rose- 
mount  858  a  5  trous  et  de  transducteurs  de  pression  pour  ihesurer  les  angles 
et  les  vitesses  des  mouvements  de  Fair  en  des  points  types  sur  un  aeronef 
de  recherche  sur  la  physique  des  nuages.  Ce  rapport  decrit  les  resultats  d’un 
essai  pousse  d’etalonnage  en  soufflerie  de  cette  nacelle  visant  a  etablir  une 
relation  entre  les  mesures  de  pression  differentielle  et  les  angles  et  vitesses 
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Symbol  Definition 


a. 
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coefficients  in  least  squares  regression  analysis 
a  constant  based  on  canister  dimensions 


pressure  coefficient 

sensitivity  factor  relating  flow  angle  to  pressure  ratio 
half  length  of  the  canister 
Mach  number 


air  pressure 

dynamic  pressure  measured  by  the  probe 
static  pressure 
total  pressure 

differential  pressure  between  the  orifices  used  for  angle  of 
attack  measurement 


differential  pressure  between  the  orifices  used  for  angle  of 
sideslip  measurement 

static  pressure  on  the  surface  of  a  sphere  at  angle  y  to  the 
free  stream  velocity  vector 
ambient  pressure 

dynamic  pressure;  wind  tunnel  dynamic  pressure 
radius  of  the  spherical  nose  of  the  canister 
free  stream  velocity;  aircraft  true  airspeed 
velocity  measured  by  the  probe 

estimated  values  of  parameters  in  regression  analysis 
longitudinal  distance  ahead  of  the  nose  of  the  canister 
longitudinal  distance  between  the  canister  nose  and  the  static 
pressure  ports  on  the  probe 

longitudinal  distance  between  the  canister  nose  and  the 

sampling  plane  of  a  PMS  2D-C  probe 

true  values  of  parameters  in  regression  analysis 

angle  of  attack 

angle  of  sideslip 

angle  between  the  free  stream  velocity  vector  and  any  point  on 
the  spherical  nose  of  the  probe 

angle  between  the  longitudinal  axis  of  the  probe  and  the 
orifices  used  for  flow  angle  pressure  measurement 
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non-dimensional  distance,  1  +  Ax/L 
air  density 
a  root-mean-square  value 

<J>  angle  of  probe  misalignment  in  roll  axi 

Subscripts 

i  incompressible 

m  measured 

o  theoretical 


1.0  INTRODUCTION 


Most  aircraft  used  in  cloud  physics  research  carry  laser 
spectrometers  externally  mounted  on  the  fuselage  or  wings  to  measure 
concentrations  of  aerosols,  cloud  droplets  and  precipitation  particles. 
Two  of  these  probes  (manufactured  by  Particle  Measuring  Systems  (PMS) 
in  Boulder,  Colorado)  can  also  be  used  to  record  particle  images  in  the 
size  ranges  25-300  ym  (2D-C)  and  200-6500  ym  (2D-P).  In  the  cloud 
physics  research  community  there  has  been  a  growing  interest  in  the 
effects  of  airflow  distortion,  caused  by  the  aircraft  and  the  pods 
housing  the  spectrometers,  on  the  concentrations  and  images  measured  by 
these  sensors.  A  workshop  convened  to  address  the  broader  topic  of 
"The  Effects  of  Airflow  Distortion  on  Aircraft  Measurements"  was  held 
at  the  National  Center  for  Atmospheric  Research  in  December,  1983 
(Baumgardner,  1984).  It  was  concluded  that  airborne  measurements  can 
be  significantly  affected  by  aerodynamic  interference  depending  on  the 
mounting  location  on  the  aircraft  and  the  parameter  being  measured,  and 
that  these  effects  are  often  subtle  and  difficult  to  quantify.  The 
scientific  community  must  give  greater  attention  to  these 
considerations  when  instrumenting  aircraft  and  using  the  resulting 
data. 

The  National  Aeronautical  Establishment  (NAE)  and  the 
Atmospheric  Environment  Service  of  Canada  (AES)  have  cooperated  for  ten 
years  to  develop  the  NAE  Twin  Otter  as  a  cloud  physics  research 
aircraft  (MacPherson  et  al,  1981,  Isaac  et  al,  1982).  The  aircraft 
contains  instrumentation  to  measure  and  record  atmospheric  state 
parameters  and  air  motion,  and  carries  four  AES-owned  PMS  probes 
mounted  on  pylons  under  the  wings  (Fig.  1).  The  effect  of  aerodynamic 
interference  on  these  probes  was  addressed  in  a  report  by  Drummond, 
1977,  followed  by  a  study  of  the  possibility  of  contamination  of  the 
measurements  by  droplets  or  ice  particles  thrown  from  the  propellers 
(Drummond,  1983).  More  recently  a  theoretical  model  for  the  flow  about 
the  Twin  Otter  wing  was  developed  that  permitted  prediction  of  the  flow 
velocities  and  angles  anywhere  ahead  of  the  wing  quarter  chord  point. 


and  in  particular  ahead  of  the  PMS  canisters  (Drummond  and  MacPherson, 
1984).  The  velocity  predictions  of  the  model  were  verified 
experimentally  by  fitting  a  standard  long  PMS  canister  with  pitot  and 
static  probes  and  pressure  transducers,  and  making  airspeed 
measurements  at  each  wing  mounting  location  over  a  wide  range  of  flight 

conditions.  The  verified  model  was  used  to  calculate  droplet 

trajectories  on  their  approach  to  the  probes  in  order  to  develop 
methods  to  correct  droplet  images  and  concentrations  for  aerodynamic 
effects  (Drummond,  1984,  Drummond  and  MacPherson,  1985). 

There  is  evidence  from  aircraft  with  wingtip  mounted  PMS 
probes  that  changing  flow  angles  due  to  formation  of  the  trailing 
vortices  may  be  causing  distortion  of  PMS  images  (Baumgardner,  1984). 
There  is  also  considerable  interest  in  verifying  the  flow  angles 
predicted  by  the  theoretical  model  for  flow  about  the  Twin  Otter  wing; 

the  original  pitot/static  equipped  canister  was  only  capable  of 

speed  measurement.  Consequently,  the  canister  has  now  been  fitted  with 
a  Rosemount  858  5-hole  probe  and  four  pressure  transducers  to  measure 
airspeed,  static  pressure  and  the  angles  of  attack  (a)  and  sideslip 
(8).  These  measurements  can  be  made  near  the  normal  sampling  plane  of 
a  2D-C  PMS  probe  or,  using  a  special  extension,  three  feet  ahead  of  the 
PMS  canister.  An  extensive  wind  tunnel  calibration  of  this  device  has 
been  performed  to  determine  the  factors  relating  the  differential 
pressure  measurements  to  flow  angles  and  velocities,  and,  in 
particular,  the  effects  of  the  canister  itself  on  these  factors.  This 
report  presents  the  results  of  the  wind  tunnel  calibration. 

It  is  expected  that  this  specially  equipped  pod  will  be  flown 
on  the  Twin  Otter  and  several  other  cloud  physics  aircraft  to  measure 
flow  velocities  and  angles  at  PMS  mounting  locations  under  wings,  at 
wingtips  and  on  fuselages.  This  report  will  serve  as  a  reference  for 
those  studies  by  describing  the  pod  and  documenting  the  calibration 
factors  necessary  for  its  use. 


A  standard  PMS  long  canister  (32£")  has  been  fitted  with  a 
Rosemount  858*  airspeed  and  flow  angle  sensor  (Fig.  2).  In  the  'short' 
configuration,  the  tip  of  the  probe  is  7.07  inches  ahead  of  the  spheri¬ 
cal  nose  of  the  canister  (Fig.  3).  This  is  about  behind  the  posi¬ 
tion  of  the  tips  of  a  2D-C  PMS  probe,  and  about  li”  ahead  of  its  laser 
sampling  plane.  This  arrangement  was  dictated  by  structural  considera¬ 
tions  and  the  desire  for  a  reasonable  spacing  between  the  canister  and 
the  static  pressure  holes  on  the  858  head.  This  configuration  will 
permit  in-flight  flow  measurements  at  a  position  representative  of  the 
actual  sampling  plane.  A  procedure  will  be  developed  below  to  correct 
airspeed  measurements  to  the  laser  sampling  plane  position. 

Flow  ahead  of  the  PMS  probes  is  of  considerable  interest 
also,  for  it  is  here  that  the  aerodynamic  forces  can  influence  particle 
trajectories  prior  to  their  arrival  at  the  sampling  plane.  Therefore, 
a  second  configuration  has  been  developed  which  employs  a  30-inch 
extension  between  the  canister  and  the  858  head  (Figs.  3  and  4).  This 
extended  configuration  will  permit  airflow  measurement  on  the  Twin 
Otter  ahead  of  the  wing  leading  edge  where  large  longitudinal  gradients 
in  flow  angle  of  attack  can  be  expected.  These  tests  will  provide  an 
additional  test  for  the  airflow  model  given  by  Drummond  and  MacPherson, 
1984. 

Aircraft  angles  of  attack  and  sideslip  are  currently  measured 
on  the  Twin  Otter  by  swiveling  vanes  mounted  on  a  noseboom  which  also 
houses  a  conventional  pitot/static  probe.  The  vanes  are  subject  to 
icing  problems  during  flight  in  supercooled  cloud.  It  is  planned  to 
replace  the  vanes  and  the  pitot/static  probe  on  the  boom  with  a 
Rosemount  858  sensor  with  a  deicing  capability.  Prior  to  this 
installation,  however,  it  was  considered  desirable  to  test  the  858 
probe  in  a  wind  tunnel  to  investigate  its  accuracy  in  sensing  airspeed 
and  static  pressure,  and  to  confirm  the  factors  relating  flow  angle  to 
the  measured  differential  pressures.  The  30-inch  extension  provides 
adequate  spacing  from  the  canister  to  allow  wind  tunnel  testing  of  the 
probe  virtually  out  of  the  aerodynamic  influence  of  the  blunt  nose  of 
the  canister,  thereby  simulating  its  installation  on  the  noseboom. 


*The  name  'Rosemount  858'  will  sometimes  be  shortened  to  '858'  in  the 
following  text. 


Figure  3  is  a  schematic  diagram  of  the  858-equipped  canister 
giving  the  principal  dimensions  for  the  short  and  extended  configura¬ 
tions.  Note  that  the  static  pressure  orifices  are  3.49  inches  behind 
the  nose  of  the  probe,  and,  in  the  short  configuration,  only  3.58 
inches  ahead  of  the  spherical  nose  of  the  pod.  Dynamic  pressure 
(airspeed)  is  measured  by  the  pressure  difference  between  the  central 
orifice  on  the  nose  of  the  probe  and  the  static  pressure  holes.  As 
will  be  shown  in  Section  3.3,  in  the  short  configuration  the  static 
pressure  will  be  influenced  by  the  blunt  nose  of  the  canister, 
resulting  in  an  underestimated  airspeed  measurement.  One  of  the  aims 
of  the  wind  tunnel  tests  was  to  determine  the  relationship  between  the 
measured  and  true  dynamic  pressures  for  later  use  on  the  aircraft. 

The  wind  tunnel  tests  were  accomplished  using  canister- 
mounted  pressure  transducers  that  were  already  available  at  the  Flight 
Research  Laboratory  and  covering  the  flight  envelope  of  the  Twin  Otter 
(145  knots  indicated  airspeed).  Table  1  lists  the  manufacturers' 
specifications  for  these  transducers.  Subsequent  to  the  tunnel  tests 
and  prior  to  flight,  it  was  decided  to  acquire  new  differential 
pressure  transducers  with  at  least  double  the  range,  so  that  a  single 
assemblage  of  transducers  could  be  flown  on  all  potential  test  aircraft 
(up  to  290  knots  IAS).  Specifications  for  these  transducers  are  also 
given  in  Table  1.  Figure  2  shows  the  pod  fitted  with  the  higher  range 
capsules  for  the  flight  configuration. 

3.0  THEORY 

3.1  Factor  Relating  Flow  Angle  to  Differential  Pressure 


TABLE  1.  SPECIFICATION  OF  PRESSURE  TRANSDUCERS  IN  POD 
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drilled  at  45  degrees  to  the  probe  axis  to  measure  flow  angles  (Fig. 

2).  For  moderate  flow  angles,  the  pressures  can  be  predicted  using 
potential  flow  theory  applied  to  the  forward  half  of  a  sphere.  In  high 
Reynold's  Number  uniform  flow  past  a  sphere,  the  static  pressure 
distribution  over  the  surface  of  a  sphere  is  given  by  Streeter,  1971: 

py  "  p„  =  *  PU2(1~  J  sin2y)  (1) 

where  y  is  the  angle  between  U  and  a  point  on  the  sphere. 

At  each  of  two  ports  on  the  sphere  at  angle  0  from  the  central  orifice, 
and  with  the  flow  at  an  angle  of  attack  a,  the  pressures  are  given  by: 


\  ~  P*  -  i  PU2(l"  |  sin2(  0+ct) ) 


p2  “  Pco  =  *  pu2(l-  f  sin2(0-cO) 


The  pressure  difference  (P  )  between  the  two  ports  is: 


— —k  =  2.  (sin2(0+a)  -  sin2(0-a)  ) 
Q  ±PU2  4 


=  -r  sin  2a  sin  20 
4 

where  Q  represents  the  dynamic  pressure. 

To  give  the  maximum  sensitivity  to  a,  0  =  ir/4  so  that: 


a  9  .  , 

~  =  4  sin  2a 


i  ,  — 1  f4  a-> 
a  =  |  sin  [j  — J 


For  small  angles  (7)  can  be  expanded  in  a  series  to  give: 


_  P  P  3  P  5 

2  r  a^  .  16  r  n\  ,  64  r 

a  "  9  J  +  2187  J  98415  '•Q  ' 


A  sensitivity  factor  kQ  is  usually  defined  so  that: 


Using  only  the  first  term  from  (8)  and  converting  a  to  degrees  gives: 


k  =  0.0785  (10) 

o 

This  is  the  value  of  kQ  normally  used  for  this  probe  for  flight  at  low 
Mach  Numbers.  The  error  in  using  only  the  first  term  in  the  expansion 
is  2%  for  Jaj<10°  and  <5%  for  Jaj  <15°. 

The  tunnel  tests  were  designed  to  confirm  the  sensitivity 
factor  kQ  by  using  the  extended  probe  configuration,  and  to  investigate 
changes  in  its  value  for  the  short  configuration  caused  by  the  presence 
of  the  blunt  nose  of  the  canister. 

3.2  Dynamic  Pressure 

A  particular  concern  in  replacing  the  conventional  pitot/ 
static  probe  on  the  Twin  Otter  with  the  5-hole  sensor  is  the  potential 
sensitivity  of  the  pressure  measurements  to  flow  angle.  Conventional 
pitot  probes  are  relatively  insensitive  to  flows  at  an  angle  to  the 
longitudinal  axis  of  the  probe,  with  errors  of  less  than  one  percent  of 
dynamic  pressure  even  at  angles  of  20°  (Ower  and  Pankhurst,  1977).  On 
the  5-hole  probe,  the  dynamic  pressure  is  measured  between  the  central 
orifice  on  the  nose  of  the  probe  (P^)  and  the  static  pressure  holes 
3.49”  behind  the  tip.  Even  at  small  off-axis  flow  angles  the 
stagnation  point  will  be  displaced  from  the  central  orifice,  and  the 
resultant  pressure  at  the  central  orifice  would  be  expected  to  be: 

P  -  p 

- —  -  1  -  ■jr  sin  o  (11) 

i  pU 

Dynamic  pressure  could  be  expected  to  be  in  error  by  up  to  one  percent 
for  flow  angles  of  only  4  degrees,  and  up  to  15  percent  at  flow  angles 
of  15  degrees.  The  second  main  objective  of  the  wind  tunnel  tests, 
therefore,  was  to  investigate  this  effect  and  determine  correction 
factors,  if  required.  An  approximate  theoretical  correction  expression 
is  developed  in  Appendix  A  which  gives: 


y  81  ^Pd  '  T  81  MPd  J  1 

y  y 


where  Pd  Is  the  measured  dynamic  presssure  P_  -  P„,  and  P  and  P„  are 
y  1  3  S’  a  B 

the  differential  pressures  for  the  angles  of  attack  and  sideslip. 


3.3  Canister  Effects 

Airspeed  is  derived  from  the  differential  pressure  measured 

between  the  central  orifice  on  the  Rosemount  probe  and  its  static 

pressure  ports.  Since  the  total  pressure  is  essentially  constant  along 

streamlines,  the  airspeed  can  be  considered  to  be  measured  at  the 

location  of  the  static  pressure  ports.  In  the  short  configuration 

these  are  only  3.58  inches  ahead  of  the  spherical  nose  of  the  canister, 

a  position  where  the  static  pressure  P  will  be  elevated  above  the  free 

b 

stream  value  P^  because  of  the  aerodynamic  influence  of  the  blunt  nose 
of  the  canister.  The  measured  airspeed  will  therefore  be  lower  than 
the  airspeed  that  would  exist  in  the  absence  of  the  canister.  The 
tunnel  tests  were  designed  to  determine  correction  factors  relating  the 
measured  dynamic  pressure  Pdy  to  the  true  values,  and  to  compare  the 
results  with  theoretical  calculations  for  this  effect. 

Three  different  theoretical  expressions  have  been  derived  to 
predict  the  effect  of  the  canister  on  the  measured  airspeed  and  dynamic 
pressure.  The  canister  was  modeled  as  a  sphere,  as  a  simple  source  and 
as  a  Rankine  body  in  a  uniform  flow  field.  The  derivations  appear  in 
Appendix  B  by  A.M.  Drummond  with  the  resulting  equations  reproduced 
below.  The  expressions  are  applicable  for  velocity  prediction  ahead  of 
the  canister  along  its  axis  of  symmetry. 

Modeling  the  canister  as  a  sphere,  the  predicted  measured-to- 
free-stream  velocity  ratio  is  given  by: 


where  R  is  the  radius  of  the  spherical  nose  of  the  canister  and  Ax  is 
the  longitudinal  distance  ahead  of  the  canister.  The  simple  source 
model  gives: 


I  =  U - i - 

u  r,  2‘M2 
(i  *— ) 


The  sphere  model  has  no  cylindrical  section  and  the  source 
model  has  infinite  length.  A  Rankine  body  has  a  finite  length  and  a 
nearly  cylindrical  centre  section  as  well  as  a  blunt  nose,  although  not 
as  blunt  as  a  sphere.  The  hemispherical  end  of  the  canister  has  a 
diameter  of  7  inches  while  the  cylindrical  portion  has  a  diameter  of 
6.5  inches  with  an  overall  length  of  32i  inches.  The  Rankine  body  can 
never  account  for  the  discontinuity  where  the  hemisphere  meets  the 
cylinder.  From  Appendix  B,  the  Rankine  Body  model  for  the  canister 
gives  the  following  velocity  ratio: 


v  0.03684  g 

"  U2-aV 


where  ?  =  1  +  —  ;  A  =  0.8989  and  L  is  the  half  length  of  the  pod 
(16.25  inches). 

All  three  models  of  the  canister  exhibit  an  axial  velocity 
defect  which  is  a  decreasing  function  of  distance  from  the  nose.  Table 
2  presents  the  velocity  ratio  as  a  function  of  Ax  for  the  three  models 


TABLE  2.  PREDICTED  VELOCITY  RATIOS 


Ax 

inches 


short  probe 


Velocity  Ratio  ■jj 


Sphere  Source  Rankine  Body 


extended  probe 


0.742 

0.879 

0.898 

0.950 

0.972 

0.982 

0.997 

0.999 


.782 

.892 

.907 

0.949 

0.968 

0.978 

0.994 

0.998 


.798 

.903 

.917 

0.956 

0.973 

0.982 

0.995 

0.998 


using  Equations  (13),  (14)  and  (15).  At  distances  greater  than  about  4 
inches  from  the  canister  nose,  the  three  models  give  almost  identical 
predictions  for  the  velocity  ratio  on  the  longitudinal  axis  ahead  of 
the  canister. 

4.0  THE  WIND  TUNNEL  TESTS 

2 

The  NAE  6x9  foot  (54.5  ft  area)  low  speed  wind  tunnel  was 
used  for  testing  the  858-equipped  pod.  The  tunnel  is  of  the  closed 
circuit  design  and  is  capable  of  steady  airspeeds  up  to  about  180  knots 
(Q  "  110  psf).  The  pod  was  mounted  in  the  centre  of  the  test  section 
upon  a  specially  manufactured  3-foot  streamlined  pylon,  which  was  in 
turn  attached  to  the  turntable  in  the  floor  of  the  tunnel.  Wiring  from 
the  transducers  was  fed  through  the  pylon  and  through  the  floor  of  the 
tunnel,  thereby  making  an  aerodynamically  clean  test  assembly.  Figures 
5  and  6  show  the  pod  in  the  tunnel  in  both  its  short  and  extended 
configurations.  The  canister  was  mounted  inverted  to  its  normal 
orientation  on  the  aircraft.  Rotation  of  the  well  calibrated  turntable 
permitted  precise  settings  of  the  angle  of  sideslip  8  of  the  probe, 
with  positive  8  corresponding  to  right  sideslip  in  the  flight 
configuration.  The  angle  of  attack  was  varied  using  interchangeable 
machined  spacing  collars  between  the  canister  and  the  pylon,  which 
permitted  nominal  angles  of  attack  of  -10,  -5,  0,  +5  and  +10  degrees. 
Angle  of  attack  was  carefully  measured  with  an  inclinometer  after  each 
collar  change. 

The  resident  data  acquisition  system  of  the  wind  tunnel  was 
used  to  record  data  from  the  pod  tests.  Output  signals  from  the  four 
pressure  transducers  were  fed  through  amplifiers  to  analog-to-digital 
(A/D)  converters,  in  which  15  bits  plus  sign  corresponded  to  ±10  volts. 
The  digitized  pressure  data  were  input  to  a  PDP-1160  processor  along 
with  measurements  of  tunnel  Q,  turntable  angle  and  tunnel  static 
pressure.  At  each  test  point  the  A/D  system  sampled  the  signal  from 
each  pressure  transducer  240  times  in  a  two  second  period.  The 
processor  computed  the  mean  and  RMS  value  for  each  pressure,  plotted 
the  four  average  pressures  versus  8  and  stored  the  output  data  on  a 


rABLE  3. 

SUMMARY  OF 

TUNNEL 

Run 

# 

Probe 

Nominal 

Q 

psf 

Ambient 

P 

OO 

mb 

Temp 

°C 

02 

short 

20  -  70 

1003 

21.7 

04 

50 

•• 

21.9 

05 

•• 

70 

- 

22.1 

07 

M 

50 

1002 

22.0 

08 

t« 

70 

•i 

22.9 

10 

•« 

50 

it 

22.6 

11 

II 

70 

22.6 

13 

•1 

50 

ii 

22.5 

14 

•  • 

70 

ii 

22.7 

16 

II 

50 

ii 

22.9 

17 

II 

70 

•• 

23.4 

19 

Long 

20  -  70 

1000 

21.9 

20 

•  i 

35 

II 

22.2 

21 

i« 

45 

•1 

22.3 

22 

•• 

55 

II 

23.9 

24 

•  I 

35 

II 

23.0 

25 

45 

11 

23.0 

26 

it 

55 

It 

23.4 

28 

it 

35 

II 

23.0 

29 

ii 

45 

II 

23.1 

30 

•  i 

55 

•1 

23.1 

32 

35 

•I 

23.0 

33 

•i 

45 

II 

23.5 

34 

•• 

55 

II 

23.5 

36 

" 

35 

•1 

23.0 

37 

•• 

45 

23.5 

38 

'■ 

55 

23.7 

8  Range 

Comments 

deg 

0 

Vary  Q 

±14 

±10 

±14 

±10 

±14 

±10 

1  ±14 

±10 

±14 

±10 

0 

Vary  Q 

±15 

±15 

±14 

±15 

±15 

±14 

±15 

±15 

±14 

±15 

±15 

±14 

±15 

±15 

±14 

disk.  These  were  subsequently  written  to  tape  and  transferred  to  the 

NRC  time-shared  IBM-3081  computer  system  for  further  processing  using 

remote  terminals  at  the  Flight  Research  Laboratory.  Figures  7  and  8 

illustrate  example  data  output  from  the  wind  tunnel  data  acquisition 

system.  The  columns  in  Fig.  7  represent  test  point,  8,  tunnel  velocity 

(fps),  P  (psf),  the  RMS  value  from  the  240  readings,  P„  and  its  RMS, 
o  8 

Pd  ,  and  the  probe  and  tunnel  static  pressure  differences  from 
ambient  (psf). 

Prior  to  installation  in  the  wind  tunnel,  the  four  pressure 
transducers  were  calibrated  at  the  Flight  Research  Laboratory  in  terms 
of  output  voltage  versus  applied  pressure.  At  the  tunnel,  the  capsules 
were  recalibrated  through  the  entire  sensor,  amplifier,  A/D  and 
processor  system  to  produce  calibration  curves  in  terms  of  pressure  in 
psf  versus  bits.  Shown  below  are  the  resulting  resolutions  for  each 
sensor  in  the  tunnel  tests  and  the  RMS  deviations  of  the  calibrated 
points  from  the  best-fitting  straight  lines. 

Resolution  RMS 

psf  per  bit  psf 

P  8.67  x  10“3  2.7  x  10~2 

°  -3  -2 

PQ  8.30  x  10  7.6  x  10 

B  -2  -1 

Pd  1.00  x  10  1.7  x  10 

y  -i 

Pg  1.64  x  10  1.3 

Table  3  presents  a  summary  of  the  test  conditions  for  each  of 
the  wind  tunnel  runs.  The  probe  in  its  short  configuration  was  first 
subjected  to  increasing  speeds  to  Q  =  70  psf  at  zero  flow  angles,  and 
then  was  tested  at  the  five  angles  of  attack  at  nominal  tunnel  dynamic 
pressures  of  50  and  70  psf.  The  extended  probe  was  tested  at  zero  a 
and  8  to  Q  *  70  psf,  followed  by  runs  with  angle  variation  at  nominal 
dynamic  pressures  of  35,  45  and  55  psf.  Table  4  summarizes  the  test 
ranges  of  speed  and  angle  which  approximate  the  flight  envelope  of  the 
Twin  Otter. 

The  test  limits  for  8  and  dynamic  pressure  were  determined  by 
structural  considerations  for  the  probe/pylon  combination.  At  high  Q, 
large  sideslip  angles  caused  vibration  of  the  3-foot  floor-mounted 


pylon.  Additional  bracing  of  the  pylon  was  accomplished  during  the 
tests,  but  the  limits  given  above  were  considered  to  be  adequate  to 
simulate  use  of  the  probe  on  the  aircraft. 


TABLE  4.  TESTED  RANGES  OF  ANGLES  AND  AIRSPEED 


Nominal  Q 
psf 

IAS 

kts 

Tested  Angle  Range  Deg. 

Short 

Extended 

a 

B 

a 

B 

20 

77 

0 

0 

0 

0 

30 

94 

0 

0 

0 

0 

35 

102 

- 

- 

±10 

±15 

40 

109 

0 

0 

0 

0 

45 

115 

- 

- 

±10 

±15 

50 

121 

±10 

±14 

0 

0 

55 

127 

- 

- 

±10 

±14 

60 

133 

0 

0 

0 

0 

70 

144 

±10 

±10 

0 

0 

5.0  DATA  PROCESSING 


A  program  was  written  to  compute  the  pressure  ratios  P^/Pd^ 


and  P./Pd  for  each  test  case  and  to  plot  these  variables  versus  the 

By 


angles  of  attack  and  sideslip.  This  program  also  employed  the  sphere 
model  for  the  canister  effects  given  in  Section  3.3  to  estimate  the 


dynamic  pressure  at  the  tip  of  the  probe  based  on  the  measured  Pd^. 


It  was  anticipated  that  the  angles  measured  by  the  probe 
might  not  be  simple  linear  functions  of  the  pressure  ratios  as 
predicted  by  Equation  (9),  particularly  in  the  case  of  the  short  probe 
where  the  aerodynamic  influence  of  the  canister  might  introduce  cross¬ 
coupling  between  the  a  and  B  measurements.  In  a  region  of  distorted 


flow,  P  ,  for  example,  might  be  a  function  of  not  only  a  but  also,  to  a 


lesser  degree,  the  sideslip  angle  B.  A  program  was  written  utilizing 


least  squares  regression  techniques  to  solve  for  the  coefficients  aQ  to 


a^  in  expressions  of  the  following  form: 


14  - 


y(i)  =  a  +  a1x1 (1)  +  a«x„(i)  ...  +  a.x.(l)  ,  i=l  to  n  (16) 
oil  L  L  j  j 

where  n  is  the  number  of  test  points,  y  is  the  'true'  value  of  a 
parameter  (such  as  tunnel  6)  and  x^  to  x.  are  the  various  estimators 
such  as  the  pressure  ratios. 

6.0  RESULTS 

6.1  Airspeed  and  Dynamic  Pressure 

As  explained  in  Section  3.3,  the  airspeed  can  be  considered 
to  be  measured  on  the  longitudinal  axis  of  the  probe  at  the  station  of 
the  static  pressure  ports.  In  the  short  configuration  this  point  is 
well  within  the  aerodynamic  influence  of  the  nose  of  the  canisLer. 

Test  Runs  02  and  19  (Table  3)  were  performed  with  the  short  and 
extended  probes  respectively,  with  a  and  6  near  zero.  The  tunnel 
dynamic  pressure  was  varied  in  10  psf  increments  from  20  to  70  psf. 
Figure  9  presents  plots  of  the  measured  Pd^  versus  tunnel  dynamic 
pressure  Q  for  both  configurations.  Clearly  the  canister  has  a  large 
effect  on  Pd^  measured  by  the  short  probe,  but  the  linearity  of  the 
plot  indicates  that  a  constant  factor  can  be  used  to  correct  Pd  to 

y 

true  dynamic  pressure.  The  least  squares  analysis  was  done  for  the 
zero  flow  angle  cases,  as  well  as  for  all  test  points  (i.e.  all  speeds, 
all  flow  angles).  The  resulting  Pd^  to  tunnel  Q  ratios  are  given 
below: 

TABLE  5.  DYNAMIC  PRESSURE  AND  VELOCITY  RATIOS 


Number  of 
test  points 

a  6 

Pdy/Q 

RMS 

psf 

v/U 

Extended-  Run  19  only 

6 

0.6  0 

0.997 

0.01 

0.998 

-  All  Runs 

251 

±10  ±15 

1.004 

0.36 

1.002 

Short  -  Run  02  only 

6 

0.2  0 

0.755 

0.10 

0.869 

-  All  Runs 

136 

±10  ±14 

0.789 

1.18 

0.888 

i '. 


Not  surprisingly,  the  least  squares  fits,  as  indicated  by  the 
RMS  values,  are  much  better  for  the  zero  flow  angle  cases.  But  for  use 
on  the  aircraft  in  flight,  there  will  in  general  be  significant  angles 
of  attack  and  sideslip,  so  the  ratios  derived  from  all  the  runs  should 
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be  used.  The  RMS  errors  of  0.36  and  1.18  psf  correspond  to  airspeed 
accuracies  of  0.4  and  1.2  knots  at  an  indicated  airspeed  of  140  knots. 
The  better  accuracy  for  the  extended  probe  is  a  reflection  of  the 
minimal  canister  effects  in  that  configuration.  Based  on  these 
results,  it  is  expected  that  use  of  the  5-hole  probe  on  the  aircraft 
noseboom  will  give  airspeeds  to  an  accuracy  of  about  0.4  knots  at  the 
speeds  and  flow  angles  experienced  by  the  Twin  Otter. 

The  last  column  in  Table  5  lists  the  experimental  velocity 
ratios  computed  from  the  square  root  of  the  pressure  ratios  Pd^/Q. 

These  are  plotted  in  Figure  10  versus  the  distance  Ax  ahead  of  the  nose 
of  the  canister,  along  with  the  predicted  ratios  from  Table  2  based  on 
modeling  the  canister  as  a  sphere,  a  source  and  a  Rankine  Body.  The 
three  models  give  very  similar  results,  and  there  is  little  aerodynamic 
effect  due  to  the  canister  for  Ax  greater  than  about  12  inches.  The 
tunnel  test  data  show  very  good  agreement  with  the  predictions.  The 
sphere  model  appears  to  provide  the  best  prediction,  and  its  ease  of 
use  also  made  it  the  preferred  model  in  studies  of  the  flow  ahead  of 
PMS  canisters  on  the  Twin  Otter  wing  (Drummond  and  MacPherson,  1984, 
1985). 

The  above  factors  relate  the  measured  velocity  to  the  true 
tunnel  velocity.  Therefore  they  can  be  used  to  correct  the  measured 
Pd  for  canister  effects  and  deduce  velocities  that  would  exist  in  the 

y 

complete  absence  of  the  measuring  pod.  In  considering  aerodynamic 
effects  on  cloud  physics  measurements,  however,  it  is  the  velocity  at 
the  laser  sampling  plane  in  the  presence  of  the  canister  that  is  of 
major  interest.  Figure  10  shows  that  at  the  sampling  plane  of  a  2D-C 
probe  the  velocity  will  be  about  0.95  of  the  velocity  that  would  exist 
in  the  absence  of  the  PMS  canister. 

The  spherical  model  (Equation  (13))  can  be  used  to  relate  the 
velocity  v^  actually  measured  by  the  858-equipped  probe,  as  derived 
from  Pd  ,  to  the  velocity  v  at  the  2D-C  sampling  plane: 


-  16  - 
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where  R  is  the  radius  of  the  canister  nose,  and  Ax  and  Ax  are  the 

P  sp 

distances  ahead  of  the  canister  nose  of  the  static  pressure  ports  on 

the  858-equipped  pod  and  the  laser  sampling  plane  on  the  2D-C  probe 

respectively.  For  the  short  probe  configuration  Ax^  =  3.58  and 

Ax  =  5.81  inches  and  the  resulting  ratio  is  1.077. 
sp 

In  summary,  when  using  the  short  probe  to  measure  velocities 
that  would  exist  at  a  2D-C  sampling  plane  with  the  PMS  pod  in  place, 
the  measured  velocities  must  be  multiplied  by  1.077.  To  estimate  the 
airspeed  that  would  exist  at  these  mounting  locations  in  the  complete 
absence  of  the  PMS  canister,  the  measured  velocity  must  be  multiplied 
by  the  inverse  of  the  experimental  ratio  given  in  Table  5,  i.e. 
1.0/0.888  =  1.126. 

As  mentioned  in  Section  3.2,  it  was  expected  that  airspeed 

(dynamic  pressure)  measured  by  the  858  probe  would  be  sensitive  to  flow 

angles.  An  expression  was  developed  (Equation  12)  utilizing  P  and  P 

a  p 

to  correct  Pd  for  potential  errors.  Considerable  evidence  from  the 
wind  tunnel  tests  shows  that  this  correction  is  generally  unnecessary, 
especially  for  the  extended  probe  which  is  representative  of  use  of  the 
858  sensor  on  the  aircraft  noseboom.  The  least  squares  regression 
analysis  using  all  test  points  was  employed  to  investigate  the  o  and  6 
effects  on  airspeed.  The  following  relationships  were  derived,  where 
Q,  the  tunnel  dynamic  pressure,  was  considered  to  be  correct. 


Q 

Q 

Q 

Q 

Q 

Q 


Extended  Probe 


0.71  +  0.996  Pd 

y 

0.71  +  0.996  Pd  +  0.0028  P  +  0.0  P„ 
y  a  g 

0.63  +  0.999  Pdy  +  0.0081  |pj  +  0.0063|pgj 
Short  Probe 
0.41  +  1.268  Pd 

y 

0.34  +  1.269  Pd  +  0.0101  P  +  0.0018  Pfl 
y  a  6 

0.51  +  1.3178  Pdy  -  0.0389  jpj  -  0.0465|Pg| 


RMS  Deviation 

psf 

0.37 

(18) 

0.36 

(19) 

0.34 

(20) 

1.18 

(21) 

1.13 

(22) 

0.51 

(23) 
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For  the  extended  probe,  the  flow  angle  coefficients  are  very 
small  and  their  inclusion  produces  no  significant  reduction  in  the  RMS 
deviations  between  the  measured  points  and  the  best-fitting  equation. 
For  the  short  probe,  the  addition  of  terms  incorporating  the  absolute 
values  of  and  improves  the  regression  fit,  with  the  RMS  error 
deviation  reduced  by  about  a  half.  Coefficients  for  the  P  and  P„ 
terras  are  5  to  7  times  larger  than  for  the  extended  case,  presumably  a 
result  of  canister  effects.  Use  of  Equation  (23)  may  be  preferred  to 
Equation  (21)  to  achieve  the  best  possible  airspeed  measurements  in  the 
short  configuration. 

An  additional  check  of  flow  angle  effects  was  conducted  by 
using  Equation  (12)  to  compute  a  corrected  Pd^  for  each  test  point. 
Figure  11  illustrates  the  corrected  and  uncorrected  dynamic  pressures 
plotted  versus  8  for  the  extended  probe  at  angles  of  attack  of  0.6,  5.7 
and  10.8  degrees.  The  tunnel  speed  was  approximately  127  knots  (Q  “  55 
psf).  It  is  apparent  that  the  uncorrected  Pd^  already  adequately 
matches  the  tunnel  Q  at  sideslip  angles  up  to  at  least  12  degrees,  even 
in  conjunction  with  angles  of  attack  up  to  10.8  degrees.  The  inclusion 
of  correction  factors  based  on  Pq  and  P^  is  unnecessary  and  would,  in 
fact,  degrade  the  Pd  measurements. 


6.2  Static  Pressure 


The  results  above  were  somewhat  unexpected,  since  total 


pressure  (and  therefore  Pd^)  was  anticipated  to  decrease  proportionally 


to  the  square  of  the  sine  of  the  off-axis  flow  angle  (Equation  (11)). 
Examination  of  plots  of  probe  static  pressure  versus  8,  such  as  shown 


in  the  lower  right  of  Figure  8,  reveals  that  P^  decreases  with  flow 


angle  in  a  manner  that  compensates  for  total  pressure  decreases,  so 


that  Pd^  is  relatively  insensitive  to  flow  angle. 


In  tunnel  testing  the  extended  probe,  static  pressure 
differences  from  ambient  up  to  4.4  psf  were  measured,  generally  in 
cases  with  high  flow  angles  at  the  higher  airspeeds.  This  could 
represent  height  errors  of  about  60  feet  at  sea  level  and  90  feet  at 
15,000  feet  when  the  858  probe  is  used  to  measure  altitude  on  the 
aircraft  noseboom.  This  performance  is  not  as  good  as  that  of  the 
conventional  static  pressure  probe  currently  used  on  the  Twin  Otter. 
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To  use  Che  858-equipped  probe  on  the  Twin  Otter  and  other 
aircraft  for  studies  of  airflow  effects  at  PMS  mounting  locations,  the 
measurement  of  a  highly  accurate  Pg  is  not  absolutely  essential  since 
this  measurement  will  likely  be  available  from  another  sensor  on  the 
aircraft.  Space  limitations  within  the  test  pod  dictated  the  use  of  a 
small  static  pressure  transducer  which  does  not  have  the  resolution  or 
the  accuracy  of  the  other  three  transducers  in  the  pod  (Section  4.0). 
Calibration  of  the  static  pressure  capsule  had  a  RMS  deviation  of 
1.3  psf. 

It  is  in  considering  the  replacement  of  the  conventional 
pitot/static  probe  on  the  noseboom  of  the  Twin  Otter  that  the  question 
of  the  accuracy  of  the  858  probe  in  sensing  static  pressure  becomes 
important.  Since  this  installation  is  simulated  by  the  extended 
configuration  in  the  tunnel  tests,  the  regression  analysis  was  used  to 
relate  the  measured  static  pressure  P  to  the  true  value.  The  best  fit 

O 

was  obtained  by  the  following  equation: 


Aps  =  Pg  -  P^  -  -2.091  +  0.1189  Pdy  -  0.0800 jP^J-  0.0684 ) P gJ  (24) 

with  an  RMS  deviation  of  0.9  psf.  The  corrected  P  could  be  deduced 
from  this  equation  to  an  accuracy  of  about  one  psf,  which  is  of  the 
same  order  as  the  accuracy  of  the  calibration  of  the  transducer.  When 
used  on  the  noseboom  of  the  aircraft,  there  are  other  position  errors 
associated  with  the  blunt  nose  of  the  aircraft  and  upwash  ahead  of  the 
wing.  Therefore,  after  installation  on  the  aircraft,  another 
expression  of  the  form  of  Equation  (24),  but  with  different 
coefficients,  must  be  derived  from  flight  test  data  using  a  trailing 
cone  or  a  tower-measured  static  pressure  as  a  reference. 

6.3  Flow  Angles 

The  theoretical  considerations  given  in  Section  3.1  indicate 

that,  for  small  flow  angles,  the  differential  pressure  measurements  Pq 

and  Pg  can  be  expected  to  be  linearly  related  to  flow  angle  in  the  form 

Pg/Q  =  kQ $ .  For  each  test  run,  that  is  at  each  combination  of  tunnel  Q 

and  angle  of  attack,  P  /Pd  was  calculated  and  plotted  versus  sideslip 

p  y 

angle  8.  The  least  squares  procedure  was  used  to  solve  for  the 


experimental  sensitivity  factor  k^  for  each  run.  Table  6  summarizes 

these  results,  and  includes  the  correlation  coefficient  (CC)  for  the 

fitted  straight  line  and  the  RMS  deviation  of  the  fit  in  terms  of 

pressure  ratio  and  angle.  Figure  12  presents  plots  of  P0/Pd  versus  8 

p  y 

for  the  extended  probe.  Each  frame  gives  results  for  a  single  angle  of 

attack  and  has  superimposed  data  points  and  fitted  straight  lines  for 

the  three  test  values  of  tunnel  Q,  nominally  35,  45  and  55  psf.  This 

figure  and  Table  6  show  that  there  is  little  variation  in  the  measured 

sensitivity  factor  k^  with  Q  or  with  a.  RMS  levels  for  the  individual 

fits  are  small,  corresponding  to  0.06  to  0.15  degrees  in  6. 

Results  for  the  short  probe  configuration  are  given  a  similar 

presentation  in  Figure  13,  with  each  frame  containing  plots  from  the 

two  test  dynamic  pressures  of  50  and  70  psf.  Here  there  is  a 

discernible  difference  in  the  slopes  for  the  two  test  velocities,  with 

Table  6  showing  that  k  increases  slightly  with  dynamic  pressure.  At 

m  ' 

Q  =  70  psf,  the  tested  8  range  was  limited  to  ±10  degrees  because  of 

the  onset  of  vibration  of  the  probe/pylon  combination  at  higher 

sideslip  angles.  The  plotted  crosses  are  well  represented  by  a  linear 

fit  with  RMS  deviations  equivalent  to  angles  of  0.07  to  0.13  degrees. 

For  Q  =  50  psf,  the  plotted  triangles  diverge  somewhat  from  the  best 

straight  line  at  angles  above  ±10  degrees.  The  resultant  RMS  levels 

fall  in  the  range  0.12  to  0.27  degrees,  still  an  acceptable  fit. 

In  Figure  14  the  test  results  for  the  extended  probe  at  all 

five  angles  of  attack  are  combined  in  single  plots.  Frames  (a),  (b) 

and  (c)  give  the  plotted  points  and  best  fitting  straight  lines  for  the 

nominal  tunnel  Q  of  35,  45  and  55  psf  respectively.  Scatter  in  the 

combined  data  is  modest,  as  indicated  by  only  slight  offsets  of  the 

overlapping  circles  representing  individual  test  points.  The  data  are 

well  fitted  by  straight  lines,  the  slopes  of  which  are  given  in  Lines 

16  to  18  in  Table  6.  Figure  14(d)  presents  data  from  all  245  test 

points,  i.e.  for  all  velocities  and  all  angles  of  attack.  The  linear 

fit  gives  k  =  0.07800  with  an  overall  RMS  level  equivalent  to  0.18 
m 

degrees  (Line  19,  Table  6). 


Figure  15  gives  a  similar  presentation  of  the  results  from 
the  short  probe  configuration.  An  average  sensitivity  factor  of  k  = 
0.1001  with  an  RMS  of  0.33  degrees  applies  when  data  from  all  130  test 
points  are  combined  (Line  32,  Table  6).  The  increased  scatter  in  the 
data  for  the  short  probe  is  likely  a  reflection  of  the  aerodynamic 
influence  of  the  canister.  Nevertheless,  the  wind  tunnel  results 
suggest  that,  for  each  probe  configuration,  the  use  of  a  single  k^  for 
all  angles  of  attack  and  all  values  of  Q  covering  the  Twin  Otter's 
speed  range  will  be  acceptable  for  flow  angle  accuracy  of  0.18  and  0.33 
degrees  for  the  extended  and  short  probes  respectively. 

The  difference  between  the  sensitivity  factors  for  the 
extended  and  short  probes  appears  to  be  due  entirely  to  the  effect  of 
the  canister  on  the  measured  dynamic  pressure  Pd  .  The  expression  for 
the  measured  sensitivity  factor  can  be  expanded  to  give: 

k  =  -  (— )  =  -  f— If-2-)  =  k  f-2-)  (25) 

m  0  LPd  >  8  lQ  jlPd  >  o  lPd  J  (  ’ 

y  y  y 

where  kQ  ■  0.0785  from  the  theoretical  development  (Equation  (10)). 
Using  the  overall  average  experimental  k  =  0.1001  for  the  short  probe, 
the  ratio  k^/k^  =  1.275.  This  agrees  to  within  one  percent  with 
dynamic  pressure  ratio  Q/Pd  of  1.0/0.789  =  1.267  determined  in  the 
wind  tunnel  tests  (Table  5),  For  the  extended  probe,  for  which  the 
dynamic  pressure  ratio  is  virtually  unity  (Table  5),  the  average  k^  for 
all  test  points  of  0.07800  is  extremely  close  to  the  predicted 
theoretical  value. 

If  the  858  sensor  head  were  to  be  mounted  at  a  different 
distance  ahead  of  the  PMS  canister,  the  results  above  show  that  the 
sensitivity  factor  could  be  reliably  predicted  utilizing  Equation  (25). 
The  pressure  ratio  Q/Pd^  could  be  calculated  by  modeling  the  canister 
as  a  sphere  and  using  the  square  of  the  resulting  velocity  ratio  from 
Equation  (13). 


TABLE  6.  LEAST  SQUARES  FITS  FOR  k 


Line  Run(s)  a 
No.  deg 


Extended 


36 

37 

38 
32 

!  33 

34 
20 
21 
22 

24 

25 

26 
28 

29 

30 

,20,24,28, 

132,26 

,21,25,39 

133.37 
,22,26,30, 

134.38 


ALL  ALL 


Short 

16 

17 

13 

14 

4 

5 

7 

8 

10 

11 

,4,7,10, 

113,16 

,5,8,11, 


Q/nom  BIAS  RMS 

psf 


79 
83 
92 
72 
.0782 
.0790 
.0767 
.0771 
.0776 
.0773 
.0783 
.0790 
.0778 
.0784 
.0792 


.0106 

.0117 

.0058 

.0098 

.0117 

.0068 

.0070 

.0076 

.0065 

.0092 

.0049 

.0123 

.0130 


.0139 


.0140 


.0156 

.0085 

.0220 

.0135 

.0268 

.0128 

.0237 

.0130 

.0119 

.0069 


CC  No.  of 

Test  Pts. 


.99997 

.99998 

.99996 

.99992 

.99990 

.99985 

.99997 

.99991 

.99985 

.99996 

.99996 

.99994 

.99996 

.99993 

.99997 


.99986 


.99985 


.99979 


.99982 


.99983 
.99991 
.99967 
.99978 
.27  .99951 

.13  .99980 

.24  .99962 

.13  .99980 

.12  .99990 

.07  .99994 

.32  .99930 

.27  .99909 


.33  .99917 


The  experimental  sensitivity  factors  k  averaged  for  each 

m 

tunnel  speed  (Table  6,  Lines  16-18  and  30-31)  show  a  slight  increase 

with  speed  over  the  tested  Mach  Number  range  of  0.15  to  0.22.  Plotted 

data  from  an  in  situ  calibration  of  an  858  probe  on  the  National  Center 

for  Atmospheric  Research  Sabreliner  showed  the  same  trend  for  Mach 

Numbers  less  than  0.35  (Fig.  4  of  Brown  et  al,  1983).  The  Prandtl- 

Glauert  Rule  (Kuethe  and  Schetzer,  1950)  indicates  that  at  Mach  Numbers 

less  than  about  0.5,  pressure  coefficients  Cp  (such  as  P„/Q)  in 

m  o 

compressible  flow  are  related  to  their  incompressible  values  Cp^  by  the 
expression: 


Over  the  range  of  Mach  Numbers  for  the  wind  tunnel  tests  of  the 
extended  probe,  a  one  percent  increase  in  the  measured  k  could  be 
expected,  which  agrees  closely  with  the  observed  values  given  in  Lines 
16-18  of  Table  6.  In  the  speed  range  of  the  Twin  Otter,  the  resulting 
variation  in  the  computed  flow  angles  would  amount  to  only  0.1  degree 
even  at  angles  of  10  degrees.  Use  of  an  average  k^  and  ignoring  Mach 
effects  would  provide  acceptable  accuracy. 

It  is  generally  assumed  that  858-type  probes  can  be  used  on 
aircraft  nosebooms  with  negligible  cross-coupling  between  the  a  and  6 
measurements  at  flow  angles  up  to  at  least  ±10  degrees  (Kaushagen, 
1949),  i.e.  that  angles  of  attack  and  sideslip  can  be  calculated 
independently.  The  tunnel  test  data  were  used  to  investigate  this 
premise  for  both  the  extended  configuration,  which  simulates  a  noseboom 
installation,  and  the  short  configuration,  where  canister  effects  might 
invalidate  this  assumption.  The  least  squares  regression  technique 
(Equation  (16))  was  used  to  relate  the  flow  angles  to  various 
combinations  of  parameters  in  an  attempt  to  determine  relationships 
that  would  provide  even  better  fits  than  the  simple  linear  expressions 
8  =*  k  ( P  / Q)  normally  used.  Results  are  presented  in  Table  7. 


Two  conclusions  can  be  drawn  from  the  equations  for  the 

extended  probe.  First,  the  inclusion  of  additional  cross-coupling  terms 

does  not  improve  the  fit  over  use  of  the  simple  linear  expression  given 

above.  The  coefficient  of  the  Pq  in  the  6  equation  is  very  small,  for 

example,  and  the  RMS  deviation  of  the  fit  is  unchanged.  Secondly,  the 

resultant  k  values  for  the  a  and  8  expressions  are  almost  identical  and 
m 

very  close  to  the  theoretical  value  of  0.0785,  which  confirms  the 
symmetry  of  the  hemispherical  end  and  pressure  ports  of  the  probe. 

In  the  case  of  the  short  probe,  the  inclusion  of  the  cross 
coefficients  appears  to  improve  the  fit  by  about  0.1  degrees  RMS.  Also 
the  k^  for  angle  of  attack  is  about  5  per  cent  less  than  that  for 
sideslip.  This  latter  effect  may  be  a  result  of  an  asymmetry  in  the 
flow  caused  by  the  pylon  supporting  the  canister.  However,  near  the 
completion  of  the  tunnel  testing  of  the  probe  in  its  short 
configuration,  it  was  discovered  that  the  858  sensor  was  improperly 
aligned  in  the  canister.  It  had  been  mounted  with  a  slight 
counter-clockwise  rotation  as  seen  from  the  front,  i.e.  the  pressure 
ports  measuring  a  were  not  perfectly  vertical.  It  was  estimated  that 
the  misalignment  was  about  2  degrees. 

TABLE  7.  REGRESSION  ANALYSIS  RESULTS  FOR  FLOW  ANGLES 


Expression 


Extended  Configuration 


8  =  0.260  +  12.8159  (Pg/Pd  )  0.0780 

8  -  0.262  +  12.8153  (Pg/Pdy)  -  0.0329  (P  /Pd  )  0.0780 

6  =  0.229  +  12.8158  (Pg/PdM  -  0.0007  Pd  y  0.0780  0.18 

a  -  -0.076  +  12.7985  (Pa/PcL)  y  0.0781  0.25 

o  -  -0.080  +  12.8052  (Pa/Pdy)  +  0.2327  (Pg/Pdy)  0.0781  0.19 

Short  Configuration 

8  -  0.142  +  9.9726  (Pg/Pd  )  0.1003  0.32 

8  -  0.144  +  9.9699  (Pg/Pdy)  -  0.2756  (Pa/Pd  )  0.1003  0.26 

a  -  0.130  +  10.4232  (Pa/Pa  )  V  0.0959  0.34 

a  =  0.134  +  10.4272  (Pa/Pdy)  +  0.3209  (Pg/Pd  )  0.0959  0.23 


RMS 

Equation 

Number 

0.0780 

0.18 

(27) 

0.0780 

0.18 

(28) 

0.0780 

0.18 

(29) 

0.0781 

0.25 

(30) 

0.0781 

0.19 

(31) 

0.1003 

0.32 

(32) 

0.1003 

0.26 

(33) 

0.0959 

0.34 

(34) 

0.0959 

0.23 

(35) 

To  examine  the  effect  of  this  misalignment,  consider  the  flow 
angles  and  8^  sensed  by  the  differential  pressure  transducers  as 
functions  of  the  true  a  and  8  and  the  rotation  angle  <J>,  such  that: 


i  pe 

8=8  cos <f>  +  a  sin$  =  ■ —  — — 
m  k  Q 

m 


(36) 


a  =  a  cos*  -  8  sin*  =  : —  — 
m  k  Q 

m 


(37) 


Solving  these  equations  for  a  and  8  gives: 


0  =  k“  ^^COS<^  ‘  k“  (Q2)Sln<f> 
m  H  m  ^ 


(38) 


01  -  (q^)sin<f> 

m  m 


(39) 


Note  that  the  signs  of  the  coefficients  of  the  cross  terms  in 
Equations  (33)  and  (35)  in  Table  7  correspond  with  those  in  the 
expressions  above.  Equating  them  gives: 

~S*n<fr  =■  -0.2756  or  =  0.3209  (40) 

m  m 

Assuming  k  =  0.100  and  solving  for  <fr,  estimates  for  <p  of  1.58  and  1.84 
m 

degrees  are  obtained,  close  to  the  estimated  misalignment  observed.  It 

is  probable  that,  had  the  misalignment  not  existed,  the  coefficients  of 

the  cross  terras  would  have  been  negligible  and  the  simple  linear 

relationships  in  Equations  (32)  and  (34)  would  provide  adequate 

accuracy.  The  misalignment  had  no  effect  on  the  determination  of  the  k 

m 

factor  from  fits  of  P^/Pd^  versus  8  since  the  cosine  of  $  is  essentially 
1.0  to  within  0.05  percent. 
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It  can  be  concluded  from  this  work  that  the  a  and  5  flow 
angles  can  be  calculated  independently  and  that  the  inclusion  of  the 
cross  terms  is  unwarranted.  For  the  extended  probe  a  sensitivity  factor 
of  0.0780  can  be  used  for  both  a  and  8,  and  the  RMS  error  levels  in 
using  the  linear  fit  to  *15  degrees  will  be  of  the  order  of  0.18 
degrees.  For  the  short  probe,  a  k^of  0.1001  appears  to  give  the  best 
results  with  an  RMS  variation  of  approximately  0.32  degrees.  It  is 
recommended  that  this  value  be  used  for  both  a  and  8,  although  the  ct  fit 
in  Table  7  indicated  a  lower  slope.  It  must  be  remembered  that  during 
the  tunnel  tests,  a  was  incremented  in  5  degree  steps  using  machined 
collars,  while  the  tunnel  turntable  provided  accurate  8  variation  in  2 
degree  steps.  The  k  factor  deduced  from  the  8  expression  should 
therefore  be  regarded  as  the  more  reliable  of  the  two. 


7.0  USE  OF  THE  PROBE  ON  AIRCRAFT 

Having  been  wind  tunnel  calibrated,  the  858-equipped  canister 
will  be  flown  on  the  NAE  Twin  Otter  and  other  cloud  physics  research 
aircraft  to  measure  airspeed  and  flow  angles  at  and  ahead  of  PMS  probe 
mounting  locations.  This  section  will  summarize  the  best  equations  and 
sensitivity  factors  to  be  used  in  airborne  tests  with  this  device  and 
present  its  operational  limitations. 

7.1  Dynamic  Pressure,  Airspeed 

For  the  extended  probe,  Pd^  =  P3  -  Pg  measured  by  the  probe 
can  be  used  to  calculate  airspeed  to  an  accuracy  of  about  0.4  knots  at 
an  indicated  airspeed  of  140  knots.  Since  canister  effects  in  this 
configuration  are  negligible,  measured  airflow  will  be  representative  of 
conditions  3  feet  ahead  of  PMS  mounting  locations  both  in  the  presence 
and  absence  of  the  PMS  probe  itself. 

Using  the  short  configuration,  the  airspeed  at  the  2D-C 
sampling  plane  will  be  1.077  times  that  derived  from  the  dynamic 
pressure  Pd  measured  by  the  probe.  To  calculate  the  airspeed  that 
would  exist  at  the  location  of  the  static  pressure  ports,  but  in  the 
complete  absence  of  the  probe  and  canister  and  their  effects,  either  of 
the  following  relationships  from  Section  6.1  can  be  used  for  the  dynamic 
pressure : 


>  "  •  *  »  •  »  m  **  •»  »  ''  *^ft  “a  *  *•  ">  •  "  *  "Tft  T«  *  «  "  >  ■  ** 


P  -  P  =  0.41  +  1.268  Pd 
T  S  y 


o  =  1.18  psf  (“1.2  knots  at  140  IAS) 


P„  -  P„  =  0.51  +  1.3178  Pd  -0.0389  P  -0.0465  P, 
T  S  y  a  ( 


a  =  0.51  psf  (“0.6  knots  at  140  IAS) 


7.2  Static  Pressure 

As  mentioned  in  Section  6.2,  better  sources  of  static  pressure 
for  altitude  calculation  will  exist  elsewhere  on  the  test  aircraft.  In 
the  short  configuration  the  pod-measured  static  pressure  will  be 
elevated  considerably  above  ambient  due  to  flow  blockage  ahead  of  the 
canister.  In  the  extended  configuration,  Section  6.2  gives  an 
expression  (Equation  (24))  with  which  to  correct  the  measured  P,,  using 
Pd^,  P^  and  P^  with  a  resulting  accuracy  of  about  1  psf. 


7.3  Flow  Angles 


The  tunnel  tests  showed  that  at  flow  angles  up  to  at  least  ±15 
degrees,  the  a  and  6  computations  can  be  handled  independently.  The 
flow  angles  can  be  treated  as  linear  functions  of  the  measured  pressure 
ratios : 


P  P 

“  =  km  (pf-")  and  8  =  km  (pf“)  (^gles  in  degrees) 


m  vPd 


For  Mach  Numbers  up  to  at  least  0.22,  the  average  values  of  the 
sensitivity  factor  k  determined  from  the  wind  tunnel  tests  provide 
acceptable  accuracy.  These  values  are  0.0780  for  the  extended  probe  and 
0.1001  for  the  short  probe,  giving  RMS  accuracies  of  approximately  0.2 
and  0.3  degrees  respectively.  At  higher  Mach  Numbers,  these  factors 
could  possibly  be  corrected  using  the  Prandtl-Glauert  Rule  by  dividing 
km  by  (1  -  M2)*. 
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7,4  Operational  Limitations 

The  differential  pressure  transducers  installed  in  the 
canister  are  limited  to  144  psf  for  flow  angle  and  288  psf  for  airspeed. 
During  the  flight  tests,  considerable  care  must  be  taken  not  to  exceed 
these  limits.  The  following  table  gives  the  flow  angle  limitations 
versus  indicated  airspeed. 

TABLE  8.  FLIGHT  LIMITATIONS 


IAS  Knots 

120 

150 

180 

210 

250 

290 

320 

Aircraft  Q  (psf) 

=  Extended  Probe  Pd 

y 

49 

76 

110 

150 

212 

285 

X 

Short  Probe  Pd  =  0.79  Q 

y 

39 

60 

87 

118 

167 

225 

275 

Flow  Angle  Limits  (deg) 

Extended  and  Short 

38 

24 

17 

12 

8 

6 

5 

The  airspeed  limits  are  320  and  290  knots  IAS  for  the  short 
and  long  probe  configurations  respectively.  The  short  probe  can 
tolerate  a  higher  airspeed  because  the  aerodynamic  influence  of  the  nose 
of  the  canister  reduces  the  sensed  dynamic  pressure.  Most  cloud  physics 
research  aircraft  make  measurement  runs  at  airspeeds  considerably  below 
these  values.  At  airspeeds  above  about  200  knots,  care  will  have  to  be 
exercised  in  limiting  the  flow  angle  excursions  during  the  airborne 
tests,  which  will  involve  Intentional  pitching  and  yawing  raanoeuvers  and 
cloud  penetrations. 

8.0  CONCLUSIONS 

The  858-equipped  canister  has  been  wind  tunnel  tested  and 
calibrated  over  the  range  of  airspeeds  and  flow  angles  representative  of 
the  flight  envelope  of  the  Twin  Otter.  The  probe  has  been  tested  in  a 


short  configuration  for  sensing  airflow  near  the  sampling  plane  of  a  PMS 
probe,  and  in  an  extended  configuration,  both  to  measure  airflow  well 
ahead  of  a  PMS  canister  and  to  simulate  the  use  of  the  5-hole  probe  on 
the  noseboom  of  the  Twin  Otter.  The  aerodynamic  influence  of  the 
canister  on  the  measured  airspeeds  and  flow  angles  has  been  documented. 

The  5-hole  probe  adequately  measures  airspeed  without  the 
necessity  of  correcting  for  off-axis  flow  at  angles  up  to  ±15  degrees. 
The  expected  reduction  in  dynamic  pressure  due  flow  angle  effects  on 
sensed  total  pressure  appears  to  be  compensated  by  matching  reductions 
in  static  pressure.  Experimentally  determined  factors  to  correct 
dynamic  pressure  and  measured  airspeed  for  canister  effects  have  been 
determined  and  agree  well  with  theoretical  predictions. 

At  flow  angles  up  to  at  least  ±15  degrees,  the  angles  of 
attack  and  sideslip  can  be  calculated  independently  without  the 
inclusion  of  cross-coupling  terms.  In  both  the  short  and  long 
configurations,  flow  angle  is  linearly  related  to  the  differential 
pressure  ratios,  with  angles  measured  to  accuracies  of  0.3  and  0.2 
degrees  respectively.  It  is  felt  that  the  relationships  determined  here 
will  still  apply  at  the  higher  speeds  on  other  test  aircraft. 

Static  pressures  measured  by  the  5-hole  probe  are  not  as 
accurate  as  those  sensed  by  a  conventional  static  probe.  In  any  case, 
when  used  on  the  noseboom  of  the  Twin  Otter,  special  flight  tests  will 
be  required  to  correct  for  position  errors  due  to  probe  effects,  the 
nose  of  the  aircraft  and  the  upwash  of  the  wing. 
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AE  TWIN  OTTER  ATMOSPHERIC  RESEARCH  AIRCRAFT 
SHOWING  THE  AES-OWNED  PMS  PROBES 


SHORT  CONFIGURATION 


FIG.  3:  SCHEMATIC  AND  DIMENSIONS  OF  THE  PROBE  IN  THE 
SHORT  AND  EXTENDED  CONFIGURATIONS 


EXTENDED  PROBE  ON  OUTBOARD  PYLON  ON  TWIN  OTTER  WING 
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FIG.  7:  EXAMPLE  PRINT  OF  WIND  TUNNEL  OUTPUT  DATA 
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FIG.  10:  PREDICTED  AND  MEASURED  VELOCITY  RATIOS 
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APPENDIX  A  -  CORRECTION  TO  DYNAMIC  PRESSURE  FOR  OFF-AXIS  FLOW  ANGLES 

When  the  stagnation  point  moves  off  the  central  pressure  port 
on  the  5-hole  probe  by  angle  a,  the  measured  pressure  would  be  expected 
to  fall  below  the  true  dynamic  pressure.  From  Equation  (1)  in  Section 


Pdy  =  P3  -  Pg  =  tPU2(1-  sin2a)  =  Q(  1—  sin2ct) 


9  .  2 


(A-l) 


The  approximation  1  -  sin  a  =  cos2a  is  accurate  to  within  2  per  cent 
at  angles  up  to  15  degrees.  Hence: 


Pd  “  Q  cos2a 

y 


(A-2) 


From  Equation  (6): 


P  =  7  Q  sin2a 
a  4 


Pg  =  Q  sin2B 


(A- 3) 


(A-4) 


Considering  first  just  the  flow  at  angle  a  (0=0),  squaring 
(A-3)  and  (A-2)  and  adding  gives: 


Pd  +  — — —  =•  Q  (sin^  2a  +  cos  2a) 
y  81 


2  16Pa  ^ 

and,  Q  “  (Pdy  +  -gj-) 


(A-5) 


(A-6) 


For  simulataneous  a  and  B  deflections,  (A-2)  would  become: 


Pd  a  Q  cos  2a  cos  20 

y 


(A-7  ) 


Squaring  and  adding  (A-7),  (A-3)  and  (A-4)  produces: 


Pd  +  — ^ —  +  ■  Q-—  “  Q  (cos  2a  cos  25  +  sin  2a  +  sin  25) 
y  oi  oi 


»  Q2 Cl  +  sin2  2a  sin2  28) 


(A-8) 


to  within  1.5  per  cent  for  a  and  8  both  10  degrees.  Therefore: 


P  2  Pi 

n  nj  Cu-  16  (  a  \  16  /'_6  'll 

Q  »  Pd  [l+  81  (pd  )  +  81  (pd  ^ 

y  y 


(A-9) 


Another  approximation,  valid  at  low  Mach  Numbers,  is  involved 
when  Q  is  represented  by  P^,  -  Pg  or  the  measured  Pd^.  From  the 
isentropic  flow  equations  in  compressible  flow  (Kuethe  and  Schetzer, 
1959): 


~  1  "  (1+ 
S 


Y  ~  1  J-  \  Y  -  1 


(A-10) 


Utilizing  the  binomial  expansion  for  the  bracketed  expression  yields: 


PT  -  Ps  -  P.ft  M2  +  l  M*  + 


(A-ll) 


But  Q  =  ipU2  =  |  Po  M2  or  P  =  ^ 


(A-12) 


and  for  air  with  y  =  7/5: 


PT  -  P  =  Q(1+t-  +  7?T  •••) 


(A-13) 


The  tunnel  tests  were  conducted  at  M  <  0.22,  so  the 
approximation  of  using  measured  P^,  -  Pg  or  Pd^  for  Q  is  accurate  to 
within  one  per  cent. 
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APPENDIX  B  -  THREE  MODELS  FOR  THE  FLOW  FIELD  AHEAD  OF  A  PMS  CANISTER 
by  A.M.  Drummond 

1.0  Introduction 

It  was  desired  to  describe  the  flow  about  the  canister  of  a 
Particle  Measuring  Systems  (PMS)  probe  using  classical  hydrodynamic 
solutions  because  the  spatial  variation  in  the  flow  velocity  in  the 
vicinity  of  the  canister  could  then  be  obtained  from  analytical 
expressions. 

The  canister  has  a  cylindrical  portion  25i  inches  long  with  a 
diameter  of  6.5  inches.  At  either  end  it  has  a  hemispherical  cap  of  3.5 
inches  radius.  The  discontinuity  at  the  joint  of  the  cap  and  the 
cylindrical  portion  can  not  be  modeled  by  simple  classical  methods. 
However,  only  the  flow  ahead  of  the  canister  along  its  axis  of  symmetry 
was  desired,  and  no  interest  was  to  be  paid  to  details  of  the  flow  aft 
of  the  canister  stagnation  point. 

Three  models  for  the  canister  were  considered  and  each  had 
certain  attractive  attributes.  The  first  treated  the  canister  as  a 
sphere  which  accurately  modeled  the  cap  but  ignored  all  effects  of  the 
canister  aft  of  the  discontinuity.  The  second  treated  the  canister  as  a 
simple  source,  which  did  not  treat  the  cap  properly  but  did  account  for 
the  finite  width  of  the  canister,  at  least  asymptotically.  The  third 
treated  the  canister  as  a  Rankine  Body  which  correctly  modeled  the 
cylindrical  portion  and  accounted  for  the  finite  length  of  the  canister 
at  the  expense  of  an  approximate  treatment  of  the  cap. 

Each  of  the  models  will  be  described  briefly  and  a  numerical 
comparison  of  the  axial  velocity  defect  ahead  of  the  cap  will  be  given 
in  the  main  body  of  the  report. 

Symbols  used  here  are  defined  within  the  Appendix  and  are  not 
necessarily  consistent  with  those  used  in  the  main  body  of  the  report. 


2.0  The  Canister  as  a  Sphere 

A  sphere  of  radius  "a"  in  a  uniform  stream  U  with  center  at 
the  origin  is  shown  in  Fig.  B-l. 


STAGNATION 


FIG.  B-1 

The  point  P  is  located  by  the  co-ordinates  r  and  0  and  qr  is 
the  velocity  component  in  the  radial  direction.  The  stream  function  is 
(Milne-Thomson,  p.  443): 


*  =  -  j  Ur2  sin2  0(l-(^)J) 


(B-l) 


The  velocity  component  is  given  by: 


1 

2  ,  a  ae 

r  sin  0 


(B-2) 


Using  <J)  from  Equation  (B-l),  the  normalized  radial  velocity  is 


given  by: 


To  consider  axial  flow,  the  angle  9  *  it.  We  define  the  point 

of  interest  to  be  Ax  ahead  of  the  stagnation  point.  The  flow  velocity 

s 

in  the  positive  x  direction  is: 

(B-4) 

(B-5) 

The  value  of  is  equal  to  1  when  Axg  =  “  and  zero  when  Axg  =  0. 


*x  ^r  ^  ^a+Ax  ^ 

s 


=  1  - 


1 


Ax  , 

-r)3 


3.0  The  Canister  as  a  Simple  Source 

A  source  placed  at  the  origin  in  a  uniform  stream  is  shown  in 
Fig.  B-2  after  Milne-Thomson,  p.  437. 


The  stream  function  is: 

i|»  ■  -  Ur^  sin^0  +  m  cos9  (B-6) 

The  value  of  m  is  chosen  such  that  the  stagnation  point  occurs 
2 

when  0  =  n  and  m  =  a  U.  Using  Equations  (B-2)  and  (B-6),  we  write  the 
expression  for  q^: 


qr  =  [cose  +  (-)  ] 


(B-7 ) 


When  0  =  it,  and  q  = 
x 


-q  ,  the  result  is: 
r 


i  ,a,.2 
q  =  1  -  (— ) 
x  r 


=  1  - 


Ax  „ 

<i+  -^)2 


(B-8) 

(B-9) 


The  value  of  "a"  is  obtained  from  the  asymptotic  value  of  y 
when  r  +  00  and  6+0,  because  y  +  R  (the  body  radius)  which  equals  2a  in 
the  limit.  Equation  (B-9)  is  rewritten  in  terms  of  R  as: 


q  =  1  - 
nx 


(B-10) 


4.  LiA.  0 

“♦  -t)! 


4.0  The  Canister  as  a  Rankine  Bod 


The  Rankine  Body  is  discussed  in  Milne-Thomson  (p.  441). 
There  is  a  source  at  the  point  +e  (Fig.  B-3)  and  a  sink  at  the  point 
-e. 

y  qr 

P  S 


r SOURCE 


STAGNATION 

POINT 

FIG.  B-3 

—e - - 

The  fineness  ratio  F  is 

defined 

as : 

„  length 

l 

(B-ll) 

diameter 

~  h 

The  stream  function  at  a  point  P  is  given  by: 


12  2 

=  —  Ur  sin  9  +  ra(cos0„-cos9.  ) 
2  Z  1 


( B— 12) 


where  9^  is  the  angle  between  the  x  axis  and  a  line  joining  P  to  -e. 
Similarly,  0^  is  the  angle  between  the  x  axis  and  a  line  joining  P  to 

e. 

At  the  point  x  =  0,  y  =  h  on  the  body,  the  polar  co-ordinates 
have  the  values : 


r  =  h 


0  -  it/ 2 


and  the  variables  in  Equation  (B-12)  can  be  written: 


ill  =  0 
0X  =  a 
02  -  Tt-a 


Hence,  from  Equation  (B-12): 


1  2 

0  =  —  Uh  -  m  2cosa 


Solving  for  m  we  get: 


4cosa 


Therefore : 


li 

2cosa 


jj—  /e2+h2 
2  e 


(B-13) 


(B-14) 


(B-15) 


(B-16) 


*,  -  %  *.  *  *  \  »„\*  *  \  O  ■  W  *  •  *  -  -  *  •  *.  •  >  •  \  *  V  -  ’ -  v.  ^  O  **  •  **  -  ■T*  "  -  n  •  *  ’  -  '  »  'i  •  *i  •  ^  O 


From  Fig.  B-3: 


2  2 

cos02  a  (rcos0-e)//r  +e  -2  ercos0 


( B— 17) 


2  2 

cos0j  =  (rcos0+e)//r  +e  +2  ercos0 


Finally,  the  stream  function  can  be  written: 


K- 


r- 


1  IT  2  .  2q  .  b  U  r 
-r  Ur  sin  0  +  -r-  (_■ 


b  U  (  (rcos 0  -  e) _  (rcos0+e) 


2  2  2  2 
/r  +e  -2  ercos0  /r  +e  +2  ercos0 


1  ,T  2  <  2a  u  h 

=  t  Ur  sin  0  +  -r  t— 
l  lie 


/  2. ,2  r(rcos0  -  e)  (rcos0  +  e) 
/e  +h  [ - r - v - 


(B-18) 


where 


J  =  /r2+e2  -  2  ercos0 


2  2 

K  *  /r  +e  +2  ercos0 


(B-19) 


We  use  the  stagnation  point  condition  to  set  the  position  of 
the  source  and  sink  inside  the  body. 

Using  the  definition  of  q^ ,  Equation  (B-18)  and  the 
simplifying  notation  of  (B-19),  we  have,  after  some  lengthy  algebra: 


5r  .  -  cos9  +  1—  /e2+h2  [(r-«  °°88)  -  (r*«  °°s6)]  ( B— 20 ) 


q.  -  sine  [1  +  T-  V’e2+li2  +  lj)] 


3  '  3- 

JJ  K 


(B-21) 


When  0  ■  0  and  r  ■  £.,  a  stagnation  point  exists  at  the  nose  of 
the  Rankine  body: 


E 


qf  =  0  =  -1  +  h  /e  +h  (-2-J-2. 

(  i  -e  ) 


(B-22) 


Equation  (B-22)  is  used  to  find  e  in  terras  of  the  fineness  ratio  F 
(Equation  (B-ll)).  We  define: 


A  =  e/£ 

and  write  Equation  (B-22)  in  terms  of  A  and  F.  With  some 
re-arrangement,  the  result  is: 


1  2  2  2 
F  (1-A  )  =  /1+A  F 


(B-23) 


(B-24) 


Solutions  to  Equation  (B-24)  for  A  as  a  function  of  F  are  best  performed 
numerically  using,  for  example,  Newton's  Method.  The  following  table 
shows  a  few  solutions. 


F 

A 

2 

.7281 

3 

.8277 

4 

.8728 

5 

.8989 

7.5 

.9330 

10 

.9499 

15 

.9666 

From  Equation  (B-20),  with  qx  =  q  ,  and  0=0  for  flow 
velocity  along  the  x  axis,  the  following  relation  can  be  written: 


9X  "  1 


-  ^-rr  TP  +  -A"  *  ■  -y^  •» 

F^  F^  ( £  -AZV 


( B-25) 


where 


Ax 

€  =  r/«  -  1  +  —■  (B- 

Note  that  the  origin  for  £  is  at  the  mid-point  of  the  body  which 
requires  that  £  be  greater  than  1. 

The  PMS  canister  has  F  equal  to  5  and  it  equal  to  16.25 

inches. 
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